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Abstract 

A new type of correlated disorder is shown to be responsible for the appear- 
ance of extended electronic states in one- dimensional aperiodic systems like 
the Thue-Morse lattice. Our analysis leads to an understanding of the under- 
lying reason for the extended states in this system, for which only numerical 
evidence is available in the literature so far. The present work also sheds light 
on the restrictive conditions under which the extended states are supported 
by this lattice. 
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The role of correlated disorder in producing extended electronic states 
in one-dimensional disordered systems has been extensively discussed in the 
literature by now The basic reason for the appearance of extended 

states in such cases has been traced to the existence of a certain type of 
short range clustering effect among the atoms, first pointed out by Dunlap 
et al. im in their study of a distribution of random dimers on a host lattice. 
It was shown that at a certain energy value the composite transfer matrix 
for a dimer offers identity contribution to the full transfer matrix, so that 
at this energy the full disordered lattice effectively behaves as an ordered 
chain. These correlations have later been shown to persist at all length 
scales in several quasiperiodic lattices like the copper mean chain and the 
period doubling lattice 0, leading to whole hierarchies of extended states. 

Instances of one-dimensional lattices are however known which offer more 
intriguing possibilities. One such example is the well-known Thue-Morse ape- 
riodic sequence, in which there is numerical evidence of extended states 
although there is no short range dimer-type positional correlations in this lat- 
tice. As we shall see in this paper, another kind of clustering effect manifests 
itself in this lattice, which cannot be analysed by the standard method known 
so far and developed in detail in Ref. |jl|] and 0. We may mention in this 
connection the recent work of Lin and Goda P] concerning a problem of hier- 
archical distribution of potentials on a lattice, where they also find an infinite 
number of extended electronic eigenstates, although there are no dimer-type 
correlations. The present work therefore provides yet another example of an 
aperiodic lattice where the standard picture of correlated disorder does not 
work. 

In a recent work on the Thue-Morse lattice, Ryu et al. ^ have provided 
numerical evidence of the existence of extended electronic eigenstates at cer- 
tain specified energies, based on the study of the trace map for this lattice. 
It has been shown that a particular value of the trace remains unchanged un- 
der a renormalisation group transformation. By using a detailed multifractal 
analysis they found that the energies associated with this trace value cor- 
respond to extended eigenstates. This work, however, does not provide any 
analytical guideline for discerning the extended character of these eigenfunc- 
tions, and the physical mechanism underlying the behaviour of these eigen- 
functions, independently of any multifractal analysis, remains to be clarified. 
The main purpose of this work is to provide this guideline. Additionally, 
we demonstrate the restrictive conditions under which this lattice supports 
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extended states. We show that extended states are supported on this lattice 
only within the framework of the on-site model; the other commonly used 
model, namely the bond-model, does not support extended states, so that 
our work invalidates the findings of Ref. P regarding the analysis of a multi- 
band model of the Thue-Morse lattice. A rigorous repudiation of the results 
of Ryu et al. has recently also been given by Bovier and Ghez . 

A Thue-Morse (TM) chain may be built up by starting with two letters 
A and B and using the inflation rule A AB and B BA repeatedly, 
where the letters A and B may be thought of as representing the atoms A 
and B. The ratio of the numbers of A and B atoms in any generation ^ is 
N a{C^) / N B^tj = 1. The electronic properties of such a lattice are generally 
modelled by the on-site version of the nearest-neighbour hamiltonian in the 
Wannier basis, the Schrodinger equation corresponding to which is given by 



n being the site index, and e„ = or e^, and the energies are in units of 
the hopping integral t which has been set equal to unity. 

The amplitudes of the wave functions at different atomic sites are deter- 
mined by using the unimodular transfer matrices of the form 



The amplitude of the wave function at the n-th site is then related to 
that at the first site by the following matrix product 



In order to unravel the correlations that are responsible for the extended 
states in a Thue-Morse chain, we first recapitulate the structural peculiarities 
of this lattice. At the very basic level, we may regard the atoms A and B to 
be the basic building blocks of the TM lattice; at the next level the pairs AB 
and BA may be regarded as the basic starting elements on which the repeated 
application of the Thue-Morse inflation rules generates the whole lattice. We 
may successively consider the pair of quadruplets ABBA and BAAB, the 
pair ABBABAAB and BAAB ABB A, the pair ABB ABAABBAAB ABBA 
and BAAB AB BAAB BAB AAB etc. as the starting blocks for generating 



{E - en)i^n = i'n+l + V'n-l 
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the Thue-Morse chain by applying the appropriate inflation rules. Evidently, 
there are no dimer or higher order atomic clustering in this lattice. It thus 
becomes necessary to analyse the products of transfer matrices corresponding 
to these strings of atoms in order to understand the conditions under which 
extended states are supported by the TM lattice. 

To handle such long products of matrices we resort to the following device. 
Let us begin by resolving the 2x2 matrices in a basis formed by the 
2x2 identity matrix I and the three Pauli matrices ax, Cy and az- We then 
have the expressions for Ma and Mb in the following form: 

Ma = oaI + Pa(^x + lA<yy + 5A(yz 

(3) 

Mb = asl + f3B<yx + Iso-y + SbcTz 

where aA(B) = 5a{b) = {E - eA{B))/2, (3a{b) = and 7a(b) = -i- 

Using Eqs.@ and the commutation rules for the matrices a^, Cy and a^, 
we may easily resolve the matrix products MaMb and MbMa to obtain 

MaMb = ail + Pia^ + licfy + 5iaz 

(4) 

MbMa = aj - Pia^ + lidy + 5iaz 

where ai = 20^05 — !, Pi = aB — aA, 7i = — ^(ayi+as) and 61 = 2aA0iB- We 
notice that the coefficients of the matrices J, ay and are the same for both 
the products MaMb and MbMa in Eqs.(^, while the coefficient of in 
these expressions appear with opposite signs. Thus if we set this coefficient, 
ie. /?! = 0, then the products MaMb and MbMa become identical, and 
the system effectively reduces to an ordered binary chain of alternating A 
and B atoms. However, it turns out that the condition /5i = can only be 
satisfied with = e^. We thus have to proceed to the next stages, and 
compare the pairs of matrix products MaMbMbMa and MbMaMaMb, the 
products MaMbMbMaMbMaMaMb and MbMaMaMbMaMbMbMa and 
so on. 

Using the basic results for the resolutions of the matrices Ma and Mb 
(Eqs.d)) and for the products MaMb and MbMa (Eqs.(|)) repeatedly, we 
can easily find the forms of the longer matrix products mentioned above. 
Interestingly, the pair of products of matrices MaMbMbMaMbMaMaMb 
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. . . and MbMaMaMbMaMbMbMa ■ ■ ■ each with 2" elements show a sur- 
prising regularity having either of the following forms 

MaMbMbMaMbMaMaMb--- = aJ + -f„ay + Sn(T, 

{n even) (5) 

MbMaMaMbMaMbMbMa--- = aj + -f'^'^y + 

and 

MaMbMbMaMbMaMaMb--- = aj + Pn^^c + InC^y + Sncr, 

{n odd). 

MbMaMaMbMaMbMbMa- ■■ = a„/ - + 7„a^ + 

(6) 

It follows that the matrix products MaMbMbMaMbMaMaMb ■ ■ ■ and 
MbMaMaMbMaMbMbMa ■ ■ ■ can never be made equal at some value of 
the energy for even values of n because their expansions given in Eqs.(|^) 
differ in two coefficients 7 and 6. On the other hand, for odd values of n, 
Eqs.(^ show that the matrix products become equal to each other if = 0. 
Thus the energy values for which = (n odd) are the ones for which 
the composite transfer matrices for of the strings of atoms ABBABAAB ■ ■ ■ 
and BAABABBA ■ ■ ■ offer identical contributions to the full transfer matrix 
for the whole chain. Since the Thue-Morse chain may be regarded as built 
out of these composite blocks, the full transfer matrix for the whole chain 
now consists of a product of identical unimodular 2x2 matrices M, each 
corresponding to a block of atoms ABBABAAB ■ ■ ■ or BAABABBA- ■ ■. 
This however does not mean that one necessarily has extended states at all 
the energies obtained by setting /3„ = 0. Several of these energies in fact 
lead to diverging, and therefore disallowed functions, and a careful analysis 
has to be made to determine precise conditions under which the string of 
such identical matrices leads to extended states, as the length of the chain 
becomes arbitrarily large. 

We first give the recursion relations connecting the coefficients a, /3, 7 
and 6 for two successive odd values of n {n > 3) which are easily obtained 
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by using Eqs.(|]) and Eqs.(|^). These are 



an = a^_2 + f3n-2 + ln-2 + + 2(7n-2 + (3a'-2 + 

+2(7'-2'^n-2 - an-2/^n-2) 

(3n = 8a„_2/3„_2(7n-2 + '^n-2) (^7) 

7n = 4(a2_2-/52_2 + 72_2 + 52_2)a„_27n-2 

Sn = 4{al_2- (3l_2 + -fl_2 + Sl_2)an-2Sn-2 

From Eqs.(0) we find that the condition /5„ = imphes either of the 
following conditions: 

an-2 = (8) 

/3n-2 = (9) 
7n-2 + C2 = 0. (10) 

First, it can be easily shown with the help of the recursion relations (0) 
that the value of = 1 in each of these cases. Among the three conditions 
Eqs.(8-10), only the condition a„_2 = automatically leads to 7^ = and 
6n = 0, while the two other conditions (9) and (10) generally yield non-zero 
values of 7^ and As will be shown presently, the first condition always 
leads to allowed values of the energies for the extended states in the TM 
lattice. The energies obtained from the two other conditions, Eqs.(9) and 
(10), generally lead to asympotically diverging functions, except in special 
circumstances to be discussed later. 

To carry out this analysis, let us express the matrix M at any level in the 
form ( suppressing the index n ) 

M=( " + ^ ^-'^ ] 
y (3 + i'y a — 6 J 

a can be identified with the quantity (l/2)Tr M. We now recall the 
well-known result for the powers of 2 x 2 unimodular matrices M 

M'" = f/„_l(x)M - Um-2{X)I (11) 
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for m > 2. Here Um{x) is the Chebyshev polynomial of the second kind, 
and X is the quantity (l/2)Tr M = a. The function Um{x) is defined as 
sin(m + 1)6/ sin 6 with cos 6' = x. In the present discussion a„ = 1 in every 
case, for which Um{l) has the limiting value (m + 1). We then have the 
following expression for M"^{x = 1): 

j^m ^ ( m{a + 6)-m + l m{/3 - i'y) 

I m{P + i'y) m(a — S) — m + 1 

Let us now choose, without any loss of generality, the amplitudes of the 
wavefunctions at the zero-th and the first site to be and 1 respectively. 
With the aid of Eq.(|l|) it then becomes simple to express the amplitudes at 
the (m + l)-th and the m-th sites as, 

i>rn+i = ni{a + 5 - 1) + 1 and ipm = m{l3 + 27) . 

We now observe that the finiteness of ipm at m ^ 00 is assured only 
when a + 5 = 1 and + ^7 = simultaneously. Now, as we have already 
seen, having a„-2 = at some stage automatically makes both 7^ and 5n 
equal to zero, in addition to making a„ = 1. Therefore, the criterion for 
the finiteness of the wavefunction at m — 00 is satisfied in this case. On 
the other hand, since neither of the conditions Eqs.(9) and (10) yield 7„ 
and 5n equal to zero, so that even though a„ = 1, the amplitudes diverge 
with increasing m. The energies obtained from these conditions may either 
not lead to proper eigenstates for the system or they may correspond to 
eigenstates with power law growth at infinity. An exception occurs when in 
condition (10) the quantities 7„_2 and (5„_2 possess common factors, so that 
energy values exist for which 7n-2 and 5n-2 become zero simultaneously. In 
this case 7^ and 5n become zero and we again have extended states at these 
energies. 

Hence we arrive at the condition q;„_2 = that gives rise to a set of 
energy eigenvalues which unambiguously correspond to extended eigenstates 
for the infinite TM lattice. Each block MaMbMbMaMbMaMaMb ■ ■ ■ or 
MbMaMaMbMaMbMbMa • • • of size 2" (n odd ) at these energy eigenval- 
ues becomes an identity matrix and therefore the total transmission across 
the chain becomes unity at these energies. For all other energy values arising 
out of the conditions (9) and (10), which of course also correspond to the 
basic condition /5„ = 0, ipm grows with m, and therefore they are not permis- 
sible. As mentioned above, the only exception occurs when 7n-2 and 5n-2 
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have common factors and therefore vanish for the same value of the energy 
at the same time. 

Although it is straighforward to obtain the algebraic expressions for the 
energy values obtained from the condition q;„_2 = 0, we calculate instead the 
energies for a specific choice of the parameter values, namely ca = — e_B = 0.5. 
We then obtain the following energy values for the first two levels : 

n = 3 : ±1.5 

n = 5 : ±2.094225204000, ±1.706503280648, 
±1.260097834749, ±0.337965671230 . 

For the sake of illustration we display in Fig.l the variation of {ipl'^ with 
the site-index n for a few selected energies. The signature of the underlying 
Thue- Morse aperiodicity is particularly apparent in Figs. 1(a), 1(c) and 1(d); 
all the wave-functions are lattice-like. The general characteristic of these 
functions is that with increasing n, amplitudes tend to cluster around groups 
of sites separated by islands where the amplitudes have very low values, a 
fact which was also pointed out by Ryu et al. 0. 

The spirit in which the present analysis has been carried out corresponds 
to the real space renormalisation group point of view. At any length scale 
corresponding to an odd value of n, Eqs.(^ define effective A and B atoms 
arranged in a TM sequence. The extended character of the eigenfunctions 
discussed in this paper owe their origin to a peculiar short range correlation 
among the atoms at every length scale, in contrast to the usual dimer-tj^e 
short range correlation discussed so far in the literature. 

The analysis so far has been made on the basis of the on-site model of 
the TM lattice. In dealing with the pure transfer model in which we have 
two different hopping strengths for the long (L) and short (S) bonds, it turns 
out that the full transfer matrix for the chain is composed of a product of 
four different transfer matrices corresponding to the LL, LS, SL and SS 
pairs. This is in contrast to the unique pair of 2 x 2 matrices that can act as 
"building blocks" for the entire TM chain in the on-site case. The four basic 
transfer matrices do not follow the TM sequence and form an intractable 
complex pattern, so that the recursive application of the TM inflation rule is 
not meaningful in this case. It therefore follows that the analysis of strings 
of transfer matrices based on Eqs.(|) and (|^) is no longer possible, and thus 
there are no extended states in the transfer model. We may point out that 
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the matrix recursion relation employed by Ryu et al. 0, M; = M/_iM;_i is 
incorrect for the transfer model, because the string of matrices do not take 
into account the atomic site flanked by a pair oi S — S bond and a pair 
of L — L bond on both sides. Here, Mi is the transfer matrix for the l-th 
generation chain, and Mi is the same matrix Mi with the roles of ti and 
ts interchanged. It can be verified easily that the above mapping does not 
generate the correct TM chain in the transfer model. Thus the analysis of 
Ref. 0] of the multiband model based on the above mapping is not valid. 
In fact, for mixed models involving variations in both site-energies and the 
hopping, we still have four different transfer matrices, and again there are no 
extended states. 

To conclude, we have demonstrated a new mechanism through which 
short range correlations lead to extended states in aperiodic one dimensional 
chains such as the TM lattice. For this, we provide a way of analysis which 
hopefully will be found useful in problems of similar kind. 
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Figure Caption 

Fig.l Plot of iV'nP versus site number n for a TM sequence with 256 sites. 
Here €a = —^b = 0.5 witli energies measured in units of hopping matrix 
element t. Figs. (a), (b), (c) and (d) correspond respectively to energies 
1.5, 0.337965671230, -1.260097834749 and -1.995507972823. 
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